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Low temperature thermal conductivity, k, of optimally-doped Bi2212 was studied before and after 
the introduction of point defects by electron irradiation. The amplitude of the linear component of 
k remains unchanged, confirming the universal nature of heat transport by zero-energy quasipar- 
ticles. The induced decrease in the absolute value of k at finite temperatures allows us to resolve 
a nonuniversal term in k due to conduction by finite-energy quasiparticles. The magnitude of this 
term provides an estimate of the quasiparticle lifetime at subkelvin temperatures. 



Fifteen years after their discovery, high-T c supercon- 
ductors continue to attract significant attention and pro- 
voke intense debate Q. One central issue is the extent 
of validity of the Fermi-liquid picture for describing the 
electronic excitations in these systems. The properties of 
the metallic state (even at optimal doping) appear to re- 
main beyond such a picture. But, are there well-defined 
quasiparticles (qp) deep in the superconducting state as 
suggested by ARPES measurements |2j? And if yes, at 
which energy scale do they break down? To answer these 
questions, low-energy excitations in the superconducting 
state are under intense scrutiny 

Low temperature thermal conductivity, k, has proven 
to be an instructive probe of such excitations. A 
non-vanishing linear term in thermal conductivity of 
optimally-doped Yf 23 for T — > was the first solid ev- 
idence for a finite density states of nodal quasiparticles 
at zero energy M. Moreover, as expected for the case 
of a d-w&ve gap the amplitude of this term was 

found to be universal; i.e. independent of impurity con- 
centration ||. Recently, Durst and Lee || showed that, 
regardless of Fermi-liquid corrections, this amplitude is 
intimately related to the fine structure of the supercon- 
ducting gap in the vicinity of the nodes. Subsequently, 
Chiao et al. [Q found a quantitative agreement between 
this gap structure as deduced from thermal conductiv- 
ity and the one directly observed by ARPES studies ||] 
in optimally-doped Bi2212. On the other hand, data 
on electronic specific heat [p|-|Tl|] and penetration depth 
|jl2| , fl3f are limited to T >2K and probe excitations within 
a higher energy interval where the density of states is a 
linear function of energy. The temperature dependence 
of electronic specific heat in Y123 is in agreement with 
theoretical expectations for the gap structure near the 
nodes gg. The same holds for the thermal variation of 



superfluid density obtained by penetration-depth stud- 
ies assuming a finite Fermi-liquid correction to charge 
transport B. Observing these remarkable quantitative 
agreements, Chiao et al. Q argued that this is a strong 
indication for the validity of Fermi-liquid treatment of 
the nodal quasiparticles. However, one important pa- 
rameter of this picture, the impurity bandwidth, 7, the 
energy scale below which the density of states becomes 
constant, is yet to be measured experimentally. 

Despite such marked accomplishment by the quasipar- 
ticle picture to treat the low-energy nodal excitations, 
numerous challenges remain. One is the absence of a 
linear term in low-temperature thermal conductivity of 
underdoped stoichiometric YBa2Cu40s H]- As there is 
no obvious reason to suppose that the structure of the 
superconducting gap in this system is radically differ- 
ent from that of the Y123 parent compound, this result 
suggests that somewhere in the underdoped regime, the 
Fermi-liquid picture might break down. Another issue is 
the effect of supercon ductiv ity on qp lifetime. Transport 

" and Bi2212 " 



studies in both Yf 23 Jl5|J]J and Bi22f 2 |3| have docu- 
mented a steep increase in the scattering time of quasi- 
particles below T c . According to recent ARPES mea- 
surements by Valla et al. jlTj on Bi2212, however, the qp 
lifetime is not affected by the onset of superconductivity. 
This surprising discrepancy remains to be explained. 

Here, we present a study of subkelvin thermal con- 
ductivity in optimally-doped Bi2Sr2CaCu20sbefore and 
after the introduction of vacant and interstitial sites by 
electron irradiation. We found that the zero-energy qp 
conductivity remains constant despite the introduction 
of pair-breaking defects, while there is a substantial de- 
crease in qp thermal conductivity, K qp , at finite temper- 
atures. This allows us to resolve the component of heat 
transport arising from finite-energy quasiparticles. 
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Two single crystals of E^S^CaCv^Os were used in 
this study with typical dimensions of 1.0 x 0.3 x 0.02 
mm 3 . Thermal conductivity was measured with a con- 
ventional two-thermometer-one-heater set-up. Point de- 
fects were introduced by exposing the sample to a 2.5 
MeV electron beam created by Van de Graaf accelera- 
tor at the Laboratoire des Solides Irradies in Palaiseau, 
France. Fig. 1 show the effect of electron irradiation 
on electrical resistivity and the corresponding low tem- 
perature thermal conductivity of the two samples. As 
seen in the two insets, irradiation leads to a decrease in 
the critical temperature and an increase in normal-state 
resistivity, p, as reported in previous electron-irradiation 
studies iQ . Both pristine samples have a transition tem- 
perature (T c =92K) and dp/dT {p = p a + bTfiQ cm, 
with b = l.5pQ cmKV 1 ) characteristic of optimally doped 
Bi2212 crystals Jl8| . In the case of sample a (p a ~60 
an irradiation flux of 6.0xl0 19 e~/cm 2 leads to a fourfold 
increase in residual resistivity (po ~240 p£l cm) and a 20 
K decrease in T c (73K). Moreover, a curvature appears 
in the temperature-dependence of resistivity which may 
suggest a doping change. In order to explore such pos- 
sibility, we have measured the room-temperature ther- 
mopower, S(300K), which is a well known indicator of 
doping levels in high-T c superconductors |^] . The deter- 
mined values of post-irradiation samples a and b (after 
the third irradiation) are -0.1 ± 0.1 and 5.0± 0.1 p V/K. 
The values of corresponding un-irradiated crystals are 
-0.8 ± 0.4 and 2.4 ± 0.1 p V/K, thus in both samples ir- 
radiation appear to have caused a slight "under-doping" . 
According to the study on S(300K) of Bi2212 @, AT C 
due to a possible variation in the doping level here are 
1.8 K and 4.5 K for samples a and 6, respectively. These 
changes are considerably smaller than AT C due to the 
pair-breaking effect from electron irradiation, 20 and 50 
K, and thus we will assume that the effect on qp trans- 
port arising from a small change in the doping level is 
negligible in comparison. It must be noted here that the 
thermal conductivity of sample b was not measured be- 
fore irradiation. Thus in the following discussions, we 
will concentrate on sample a and will use the data on the 
sample b only as supplementary evidence. 

The main panels of Fig. 1 present the low temperature 
thermal conductivity data. The persistence of a size- 
able phonon contribution, n p h, to the total heat trans- 
port complicates the determination of qp conductivity. 
As T — > 0, K p h, is expected to enter the boundary scat- 
tering (or ballistic) regime, where phonon mean-free-path 
is limited by the dimensions of the crystal and n p h be- 
comes proportional to T 3 . Above this regime, K p h should 
increase more slowly due to the decrease in phonon mean- 
free-path. Therefore, by plotting n/T against T 2 , one 
can resolve the zero- intercept, k^q/T, corresponding to 
the residual linear component of k that is a signature 
of residual qp conductivity. As seen in fig. lb, there 
is a finite kqo/T term in the thermal conductivity of 
sample a both before and after irradiation. The mag- 
nitude of this term shows little or no change in spite of 



the drastic increase in the number of defects. Note that 
the original contacts were kept throughout the succes- 
sive measurements for each crystal, the sizeable uncer- 
tainty in k geometric factor does not hamper the effect 
of disorder. If we assume a linear extrapolation of the 
normal-state resistivity, the scattering rate in sample a 
becomes four times larger after irradiation. If we further 
take the enhancement in the qp lifetime below T c into 
account, the induced increase in impurity scattering rate 
becomes even larger. The amplitude of kqo/T found in 
sample a (0.16 ± 0.03 mW/K 2 cm) is similar to that of 
sample b (0.13 mW/K 2 cm) as well as to those reported 
in previous studies (0.14-0.15 mW/K 2 cm) |2^,g|. Theory 
predicts an increase in kqq/T at sufficiently high levels 
of impurity concentration due to the impurity-induced 
change in T c p2]| . Therefore, considering that the sam- 
ple b after the third irradiation is nearly twenty times 
dirtier than the pristine sample a (judging from their re- 
spective Po values), the "universal" character found in 
the values of Koo/T among all samples here is quite sur- 
prising. A quantitative comparison to the fine structure 
of the superconducting energy gap around the nodes can 
then be performed by inserting the value of kqo/T into 
the following equation || : 

/rr k B 2 v F n 

koo T = — -. (1) 

in V2 d 

where vf and i>2 are qp Fermi velocity and gap velocity at 
the nodes in the gap. One then obtains ^ = 21±3 which 
is close to the ratio obtained by ARPES measurements 
||. The same conclusion has been previously drawn by 
Chiao et al. {§. 

Before discussing the finite temperature qp conductiv- 
ity further, let us focus on phonon conductivity. For the 
sample a, if one can assume that the system enters the 
ballistic regime for T < 0.25 K, the extracted cubic term 
for phonon conductivity is n p h/T 3 ~ 2.6mW/K cm. 
This can be compared to the theoretically expected 
value using k p h = l/3/3<v p h>l p h,T 3 , where /3 is the 
phonon specific heat coefficient and v p h is the aver- 
age sound velocity. Taking the transverse dimensions 
of our crystal along with the available values for (3 — 
0.0095mJ/K 4 cm- 3 p§ and v ph = 3600ms" 1 |§, n ph /T 3 
should be 5.9 mW/ IC^cm, approximately twice the value 
determined from our data. One possible source for this 
discrepancy is that our measurement stops above the on- 
set temperature of phonon ballistic conductivity. In the 
study reported by Chiao et al, the ballistic regime ap- 
pears only below 130 mK. In this downward cur- 
vature of the k/T curve would lead to a slight decrease 
in the finite intercept and the estimated koo would be- 
come 0.14 mW/K 2 cm (see Fig. 2a) but remains within 
our experimental uncertainty. 

Another important aspect of Fig. lb is a sizable 
irradiation- induced decrease in k(T) for the entire tem- 
perature range (0.13 K <T < 0.9K). We begin by not- 
ing that the typical wavelength of acoustic phonons may 
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be estimated to vary as X p h = hvph/ksT — 173nm/K. 
Thus, at subkelvin temperatures, the spatial extension 
of lattice vibrations is more than two orders of magni- 
tude larger than the size of introduced point defects and 
thus Kph should not be affected by irradiation. For this 
reason, the observed decrease must be exclusively due 
to an increase in the qp-defect scattering rate. However, 
the qp conductivity is expected to become independent 
of impurity scattering rate for T<7, the impurity band- 
width. Such condition is finally met in sample b after 
the second and the third irradiation where k is virtually 
unchanged (see the lower panel of Fig.l). In relatively 
cleaner samples, the marked decrease in k infers that 7 
must lie below our range of measurements. 

Fig. 2 depicts the change in n(T)/T of sample a be- 
fore and after irradiation. Ak/T is linear in T, revealing 
a quadratic temperature dependence. The most plau- 
sible explanation for this result is to concede that qp 
conductivity of the pristine sample includes a quadratic 
term which is heavily diminished by irradiation. A fit 
to the data of Fig. 2 for the whole temperature range 
yields An = aT + bT 2 with a =0.005mW/K 2 cm and 
b =0.19mW/K 3 cm. The small size of the linear intercept 
indicates again that irradiation has left the linear term of 
qp conductivity intact. We now compare the amplitude 
of the term b with what is expected from the finite-energy 
qp contribution. For energies exceeding 7, the density of 
states varies linearly with energy I^Jij: 
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This leads to a quadratic temperature dependence of spe- 
cific heat U 
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where £(3) ~ 1.20 is a numerical factor. Thermal conduc- 
tivity of these excitations can be estimated via kinetic 
theory; K e =|C e u| r e . Here, r e is the electronic scatter- 
ing time. The temperature dependence of r e governs the 
temperature dependence of electronic thermal conduc- 
tivity. The relative weight of this term to the universal 
linear term can be estimated to be: 



Ke_ = 18C(3)fc B T e T 



(4) 



Note that while koq/T is universal, n e is not: its 
magnitude decreases with increasing disorder. The 
size of the quadratic term coefficient appearing in 
Ak(T) (b =0.19mW/K 3 cm) corresponds to an electronic 
scattering time of 1.3 ps. On the other hand, qp life- 
time can be estimated from p — — and w„ = 

(47rne 2 /m*) 1 / 2 , where uj p is the Drude plasma frequency. 
Using uj p = l.leV & p(100K) = 200/ificm and p = 
one finds t^IOOK) ~ 0.05 ps and T e (0K) ~ 0.2 
ps, provided that the scattering rate is also linear in tem- 
perature for T < T c . Our results suggest a significant 



increase in qp lifetime, and thus, a substantial reduction 
in scattering rate below T c , in agreement with microwave 
conductivity data It is instructive to compare the 

magnitude of r e in Bi2212 and Y123. In the latter sys- 
tem, T e of high-quality crystals is of the order of 7 ps 
fl5f and is reported to approach 20 ps ||(J in BaZrOa- 
grown crystals. For Bi2212, recent study on the complex 
conductivity by Corson et al. showed the qp lifetime to 
approach ~1 ps |27j] at low temperatures, similar to the 
value calculated independently from our measurements. 
Further, assuming the scattering rate to remain constant 
p6| , one can estimate the size of K e upward in tempera- 
ture. The deduced value at 5K (~4.8 mW/Kcm) can be 
compared via the Wiedemann- Franz law to the available 
data on charge conductivity which is limited to T >5K 
fU§. This yields ^f(5K) ~ 21.3 nW/fi K, very close to 
the Sommerfeld value (Lo=24.5 nW/fi K). Thus, our in- 
terpretation appears to be consistent with what is known 
from charge conductivity. 

In spite of this apparent consistency, this analysis can- 
not accommodate the broader theoretical picture of qp 
transport in d-wave superconductors elaborated during 
the recent years. First, a linear N(E) implies an energy- 
dependent r e . In the unitary limit, for example, this leads 
to a cubic (instead of a quadratic) behavior for finite- 
energy qp transport for T > 7, @. Second, the size of 7 
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strongly depends on the impurity density 
as the scattering phase shift 8. In the unitary (6 = 7r/2) 
limit, 7 increases substantially with increasing scattering 
rate, T; 7 ~ (rA ) 1/2 . In the Born (S = 0) limit, its 
enhancement with T is largely attenuated for T <C T c ; 
7 ~ Aqcxp(— ^t-). Here V = and Ao is the magni- 
tude of the superconducting gap || . Now, with r e ~ 1 
ps (which would yield T ~ 3K) and A - 40meV 
7 may be estimated to be 35K in the unitary limit and 
virtually zero in the Born limit. Thus, the size of 7 found 
in this study is in sharp contrast with what is expected 
in the unitary limit. A slight deviation from the unitary 
limit can produce a linear k(T)/T above its universal 
value in a limited temperature range [ p9[ , however, pre- 
liminary investigations reveals that a reasonable phase 
shift only cannot account for the magnitude of the ex- 
cess conductivity observed here |3(| . Clearly, our findings 
constitute a challenge to existing theory. 

In summary, we have studied the effect of electron ir- 
radiation on the low temperature thermal conductivity 
of optimally doped Bi2Sr2CaCu20s- The quasiparticle 
contribution to heat transport was found to contain two 
distinct components, a linear term associated with zero- 
energy quasiparticles and a quadratic term originating 
from finite-energy quasiparticles. The linear term re- 
mains insensitive to the number of defects. The mag- 
nitude of the quadratic term allowed us to make a new 
independent estimation of quasiparticle scattering time 
at low temperatures, implying a significant increase in 
the quasiparticle transport lifetime below T c as indicated 
by other probes. 
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FIG. 1. Upper panel: Low-temperature thermal conduc- 
tivity, k/T, of sample a is plotted as a function of T 2 . The 
thin lines are guides to the eye. The thick line represents 
the expected asymptotic lattice conductivity at the ballistic 
regime when the phonon mean-free-path attains the average 
sample size (see text). The inset shows the resistivity data on 
the same sample. Lower panel: Same for sample b. Note that 
the theral conductivity of this sample in the pristine state was 
not measured. 



however, that according to these authors, the onset of 
the superconductivity does not affect the scattering rate. 
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FIG. 2. Upper panel: The change in the thermal conduc- 
tivity of the sample before and after electron irradiation di- 
vided by temperature and plotted as a function of tempera- 
ture. The straight line is a fit to the data revealing a neg- 
ligible intercept and a quadratic variation of Ak(T). Lower 
panel: Thermal conductivity of the sample before and after 
irradiation together with a sketch of electronic and lattice 
components of thermal conductivity in the pristine sample. 
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